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Abstract. We propose a model describing liquid-solid phase coexistence in mixed lipid membranes by 
including explicitly the occurrence of a rippled phase. For a single component membrane, we employ a 
previous model in which the membrane thickness is used as an order parameter. As function of temperature, 
this model properly accounts for the phase behavior of the three possible membrane phases: solid, liquid 
and the rippled phase. Our primary aim is to explore extensions of this model to binary lipid mixtures 
by considering the composition dependence of important model parameters. The obtained phase diagrams 
show various liquid, solid and rippled phase coexistence regions, and are in quantitative agreement with 
the experimental ones for some specific lipid mixtures. 



PACS. 87.16.Dg Membranes, bilayers, and vesicles. 

1 Introduction 

In recent years, domain formation in biomembranes and 
artificial membranes has attracted great attention in con- 
nection to the so-called "raft formation" in biological cell 
membranes [Tj. It is beheved that the appearance of such 
domains (rafts) in membranes plays an important role for 
various cell functions pl[3]. Beside their interest in biocel- 
lular processes, raft formation has a fundamental physi- 
cal interest, because it offers an special example of two- 
dimensional phase separations coupled to internal and ex- 
ternal membrane degrees of freedoms Using fluores- 
cence microscopy or X-ray diffraction techniques, the lat- 
eral phase separation between different Hquid phases, or 
between the solid and liquid phases has been observed for 
ternary-component vesicles consisting of saturated lipid, 
unsaturated lipid, and cholesterol [HinilTlIHlIS] . Although 
the biological signiflcance of rafts is not yet fully under- 
stood, it serves as our primary motivation to explore, from 
a physical point of view, phase transitions, domains and 
phase coexistence in multi-component lipid membranes. 

Lateral phase separation in membranes occurs even 
for the simpler case of binary lipid mixtures without any 
added cholesterol. The coexistence is between liquid-like 
(La) and solid-like (L/j' called "gel") phases, and it can 
be visualized using techniques such as two-photon fluo- 
rescence microscopy [lO]. For example, for DPPC/DPPE 
Hpid mixture, the solid domains exhibit morphologies such 
as hexagonal, dumbbell or dendritic shapes [10], while for 
DPPC/DLPC lipid vesicle the soHd domains appear to be 
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stripe-like and are monodispersed in their domain width 
[llj . Shapes of solid domains found on spherical fluid vesi- 
cles has been also addressed theoretically [12] • By assum- 
ing the additivity of stretching and line energies, the model 
predicts a phase diagram including cap, ring and ribbon 
phase domains on a spherical surface. 

Although visualization of domains embedded in fluid 
membranes has become possible only recently, the phase 
diagram of various binary lipid mixtures in bilayer mem- 
branes has been known for some time [ISjTiJ. In a pre- 
vious model investigated by some of us [15,16j, a cou- 
pling between composition and internal membrane struc- 
ture was proposed and the resulting phase diagram was 
calculated. It includes a coexistence region between La 
and L(3> phases in agreement with experiments on speciflc 
lipid mixtures. However, for some lipids, such as DMPC or 
DPPC, another distinct soHd-Hke phase is known to occur 
and is called the "rippled" or Ppi phase. The unique feature 
of this phase is that the membrane shape is spatially mod- 
ulated while the lipid hydrocarbon tails are ordered. The 
speciflc phase diagrams of DMPC/DPPS, DPPC/DPPS 
[T7] . or DPPC/DPPE [18| mixtures have been explored 
in experiments, and show the three phases {La, Lpi, Ppf) 
together with the coexistence regions between them. How- 
ever, these phase diagrams have not yet been considered 
theoretically. 

The main objective of the present paper is to provide 
a simple model that describes the liquid-solid coexistence 
in binary Hpid mixtures while considering the possibility 
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of an intermediate rippled phase. Our starting point is a 
model introduced by Goldstein and Leibler (GL) that ac- 
counts for the succession of phase transitions as applied 
to single-component lipid membrane [T9l[20] . In the GL 
model, the membrane thickness is used as a scalar order 
parameter. We combined some of the ideas presented in 
the GL model with our own previous model for binary 
(or even ternary) lipid mixtures |15lil6j . For binary lipid 
mixtures, the model parameters are taken to be depen- 
dent on the relative lipid composition, and resulting in 
several types of phase diagrams that are in quantitative 
agreement with experiments on specific lipid mixtures. 

In the next section, we first review the GL model de- 
scribing the phase transitions in a single-component lipid 
membrane. The corresponding mean-field phase diagram 
is obtained in terms of the temperature and membrane 
elastic constants. In Sec. [HI we propose an extension of 
the GL model for a two-component lipid membrane, and 
calculate the binary phase diagrams. Some discussions and 
comparison to other models are discussed in Sec. [H 



2 Single-component lipid system 

In a single-component lipid bilayer, one typically observes 
a discontinuous first-order phase transition from the La 
to the phase as the temperature is decreased. This 
is called the ^'main transition" and is associated with or- 
dering of the hydrocarbon tails of the lipid molecules. In 
the La phase, the hydrocarbon chains are disordered (and 
the phase is liquid-like), while they order by stretching 
and tilting in the L^' (solid-like) phase. Several theoreti- 
cal models have been proposed in the past to describe the 
main transition in isolated membranes |2ip22j. 

In other single-component lipid systems, however, a 
rippled {PfS') phase appears between the La and Lp' phases 
for stacked bilayer membranes [23jl24jl25| . This phase is 
more peculiar and many studies have been devoted to 
understand better its modulation wavelength and ampli- 
tude [26]. For a stack of lipid bilayers in water (a lamel- 
lar phase), the phase diagram as a function of relative 
humidity and temperature was reported [271128] . Several 
attempts have been made to describe the rippled phase 
theoretically such as those based on a molecular level de- 
scription [291130] . or Monte Carlo simulations [31]. In con- 
tinuum theories of the rippled phase, different quantities 
have been suggested for the order parameter. Examples 
are the membrane thickness [19|l32j or the configuration 
of the hydrocarbon chains (ratio of the trans bonds in the 
chains) [33]. In both cases the order parameter is scalar. 
More recently, a Ginzburg-Landau theory was proposed by 
Chen, Lubensky, and MacKintosh (CLM) who employed a 
vectorial order parameter representing the tilt of the lipid 
molecules [34|l33]. 

We further discuss the CLM model in Sec. |4] below, 
but our own starting point is based on a model proposed 
by Goldstein and Leibler (GL) for single-component lipid 
membranes |19p20j. To describe the main transition in- 
volving chain ordering and stiffening, GL introduced a 



dimensionless scalar order parameter of the membrane, 
m{r): 

6{r) - So . 

m{r) ^ . (1) 

This parameter depends on the actual membrane thick- 
ness (5(r), and on the constant membrane thickness i5o of 
the La phase. The two-dimensional lateral position within 
the bilayer plane is denoted by r. Notice that m encapsu- 
lates changes that may occur in several degrees of freedom, 
including, for example, the conformations of the hydro- 
carbon chains, molecular tilt and positional ordering. In 
the present paper, the bilayer nature of the membrane is 
not taken into account and the bilayer thickness is simply 
taken as a sum of the two monolayer thicknesses. The role 
of the bilayer structure in the formation of rippled phases 
was explicitly considered in previous models [36|137|. 

For simplicity, we assume that the rippled phase is 
spatially modulated only in the a;-direction, and ignore 
two-dimensional rippled phases such as the square lattice 
phase [34ll35] or the hexagonal phase |i38.:39j . For an iso- 
lated lipid bilayer membrane, the stretching free energy 
per lipid molecule is [19] 
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A very similar free energy was proposed by Marder et 
al. [32J. The first three terms are the Landau expansion 
in powers of the order parameter m. Only 02, the sec- 
ond order term coefficient, has an explicit temperature 
dependence: 02 = a'2(T — T*), with T* being a reference 
temperature. (It actually is the critical temperature in 
the absence of the cubic term.) Because the Lp' La 
phase transition is known to be first-order, 03 is taken to 
be negative, whereas 04 is always positive to ensure sta- 
bility of this free energy expansion. The next two terms 
are related to the lowest order gradients of m. (Note that 
the x-coordinate is rescaled in units of a molecular length 
scale £ so that x x/i is dimensionless.) These gradi- 
ent terms represent elastic out-of-plane undulations of the 
membrane and their coefficients are the elastic constants 
C and D, respectively. The coefficient C can be either 
positive or negative, but D is always positive in order to 
ensure stability of the expansion. The physical origin of a 
negative C value can be related to the coupling between 
the conformation of the chains and the curvature of the 
lipid/water interface, or to interactions between the polar 
head-groups of the lipids and water |19|. This point will be 
further discussed later in Sec. HI In the absence of spatial 
gradient terms in Eq. ^ and if m{x) is constant in space, 
the main transition temperature Tm is calculated from the 
conditions fst = d/st/dm = 0. It is related to T* by 
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With nonzero 03, a first-order melting transition at T„ 
preempts the critical point at T*. 
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Fig. 1. Mean-field phase diagram of a single-component lipid 
membrane as a function of (negative) elastic constant C and 
temperature T (in degrees Celsius). The liquid (ia), solid {L(}i) 
and rippled {Ppi) phases are separated from one another by 
first-order phase transition lines which meet at a triple point 
(marked by a full circle): Ctr = -3.33 x 10"" J, = 38 °C. 
The parameters used to calculate this phase diagram are given 
in the text. 

When C > 0, the equihbrium phase is always homoge- 
neous in space; either as a Lfji or La phase. On the other 
hand, for C < 0, the rippled phase with a characteristic 
wave vector q* = {—0/20)^/'^ (in unit of may dom- 
inate over the two other phases. (It should be noted that 
C < is only a necessary condition to have the rippled 
phase as will be discussed below.) To obtain the mean-field 
phase diagram, we employ the single-mode approximation 
in which the rippled phase is described by the most domi- 
nant wavevector q* , only. This approximation is valid near 
the phase transition point where the segregation tendency 
between the two lipids is weak enough. The order param- 
eter m can be written as |32j 

m = mo + TOi cos((7*a;), (4) 

where mo = (m) is the spatial average of m, and mi is 
the amplitude of the single q*-moAe of the rippled phase 
taken to lie arbitrarily in the x-direction. 

By substituting Eq. ^ into Eq. ^ and taking the 
spatial average over one period, we obtain 

(f^t) — fTi? H — flo I mn H — m? I 

\jst/ -^g^ 1 2 V 2 V 

1 3 3 2\ 

+ 3^3 I mo -f -mom^ 1 

-|~-a4 (niQ + Smlmf + ^"^i ] ■ (5) 



The two homogeneous Hquid and solid phases, Lq and L/ji, 
respectively, are characterized by mi = 0; while mo = 
in the disordered liquid La phase, and is non-zero in the 
soHd Lp> phase. From the three phases only the P/3' phase 
is characterized by a modulating amplitude, mi 7^ 0. 

The mean-field phase diagram is obtained by minimiz- 
ing Eq. 1(5]) with respect to both mo and mi, and compar- 
ing the relative stability of the three relevant phases. In 
Fig. [U we present an example of such a phase diagram, 
as function of the (negative) effective elastic constant C 
and temperature T plotted in degrees Celsius. The param- 
eter C can be a function of relative humidity (fraction of 
water content in the lamellar phase) as discussed below. 
The phase diagram is calculated for a choice of system 
parameters which reproduce the main transition temper- 
ature of DMPC at = 38 °C: 4 = 2.4 x IQ-^ijK-i, 
ag = -1.1 X 10-i«J, a4 = 2.2 x IQ-^^^J, T* = -13 °C 
(260 K), and L» = 2.0 x lO-^^J. These values (except 
D > whose value is less important) are taken from 
Ref. |2Q]. By inserting the above Landau coefficients into 
Eq. ^ , it can be readily checked that the main transition 
temperature of DMPC = 38 °C (Fig. [J) is obtained. 

In Fig. [J, three first-order phase transition fines sep- 
arate the three phases and meet at a triple point, Ctr = 
— 3.33x 10~^^ J. Due to the first-order nature of the L/j' — > 
La phase transition, it can be argued, on general grounds, 
that the triple point Ctr has to be located at non-zero 
values of C. For small magnitude of C, |C| < |Ctr|, the 
horizontal first-order phase transition line occurs at T = 
Tm = 38°C and is C independent. An increase in temper- 
ature (as long as |C| < |Ctr|) will melt the sofid phase Lp' 
directly into the liquid phase La at Tm. But for negative 
and large enough magnitude of C, |C| > |Ctr|, any in- 
crease in temperature will cause a sequence of phase tran- 
sitions: first the solid phase Lpi melts into the P/s' phase, 
and only then, upon further increase of the temperature, 
the P/3' phase will make a phase transition into the liq- 
uid phase La- It is also apparent from the figure that the 
region of the P^' phase expands on the expense of the uni- 
form liquid and solid phases as \C\ increases. This means 
that the phase transitions L/si Ppi and P/3' La (at 
constant C) occur at larger temperature deviations from 
Tm as |C| increases. 

The phase diagram of Fig. [1] essentially reproduces 
all the experimental facts observed for single-component 
lipid membranes. It qualitatively agrees with the observed 
phase diagram of a lamellar phase of DMPC, when the 
negative elastic constant C is taken to be proportional to 
the relative humidity ^27u28 j. Although the exact depen- 
dence of the elastic constant C on humidity is not known, 
the former can possibly be reduced by increasing hydra- 
tion. For another lipid, DPPC, the Pp' phase is found 
experimentally to occur at temperatures around 37 °C (in 
Ref. [18] the precise temperatures for the L^i ^ P^i and 
Ppi — > La phase transitions have not been fully reported) , 
while the rippled phase is not observed for a stack DPPE 
lipid bilayers. Such a difference in the Hpid phase behavior 
can be qualitatively attributed to different values of the 
elastic constant C. 
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Fig. 2. Calculated mean-field phase diagram of a binary 
lipid mixture as a function of their relative composition 4> 
and temperature T. The parameter values are chosen to fit 
DPPC/DPPE mixtures {4> = Q represents pure DPPC and is 
denoted by a subscript "B"): Ca = -2.0 x 10"^^ J, Cb = -5.2 x 
10"^^ J, TX = 9 °C, Tb = -15 °C, and J = 1.45 x 10"^" J. All 
other parameter values and definitions of the different phases 
are the same as in Fig. [U The critical point is indicated by a 
full circle and occurs at Tc — 21.3 °C. The horizontal dashed 
line indicates the three-phase coexistence at the triple point, 
Ttr = 42.0 °C. 



3 Two-component lipid mixtures 



In this section, we extend the above single-component 
model to membranes consisting of a binary lipid mixture. 
We employ a similar approach as was used in our previ- 
ous study [TSlfTB] . and consider the coupHng of the melt- 
ing phase transition with a lateral phase separation in the 
mixed membrane. 

A binary mixed membrane is modeled as an incom- 
pressible A/B mixture of 4> mole fraction of lipid A and 
{\ — 4>) fraction of lipid B. For simplicity, we assume the 
same area per molecule for both species and ignore any 
Hpid exchange with the surrounding solvent. In general, 
the two lipids will have different main transition tempera- 
tures originating from different molecular parameters such 
as chain length, degree of saturation and hydrophilic head 
group. 

The total free energy per Hpid, /tot = /st + /mix, com- 
prises: (i) the chain stretching free energy /st given by 
Eq. |[2]), and (ii) the free energy of mixing, /mix- The lat- 
ter energy per Hpid molecule is the sum of the entropy of 
mixing and enthalpy. It can be written within the Bragg- 
Williams (mean-field) approximation as 

/mix = k^T[c^ \ogc^+ (1 - 0) log(l - 0)] -f i J0(1 - 0), (6) 



Fig. 3. Calculated mean-field phase diagram of a binary lipid 
mixture as a function of their relative composition and tem- 
perature T. The parameters are the same as in Fig. [2] except 
J = 1.85 X 10"^" J. The critical point is located at 62.1 °C, 
and three triple points occur at 35.6, 41.3, and 56.3 °C, respec- 
tively. This choice of model parameters is in agreement with 
experimental phase diagram of DEPC/DPPE mixture where 
= represents pure DEPC. 



where fee is the Boltzmann constant, and J > is an 
attractive interaction parameter between the lipids that 
enhances lipid-Hpid demixing. 

For a binary mixture the free energy /st is assumed 
to have the same functional dependence on the effective 
elastic constant C and the reference temperature T* . Al- 
though these two parameters depend on the lipid compo- 
sition (j), the precise dependence cannot be calculated from 
such a phenomenological approach. Alternatively, we pro- 
ceed by further assuming the simplest linear interpolation 
between the two pure lipid limiting values: 



= 0Ca + (1-(/')Cb, 

T*(0) = 0Ti + (i-0)r^, 



(7) 

(8) 



where Ca and Cb are the elastic constants, and and 
Tg are the reference temperatures of Hpid A and B, re- 
spectively. The other elastic constant Z3, as well as 03 and 
a4 in Eq. ([2]) are assumed to be the same for the two lipids 
because the important properties of the lipids are mostly 
reflected in the choice of their C and T* parameter val- 
ues. As for all other model parameters, we use hereafter 
the same values listed in the previous section. 

The total free energy /tot = /st + /mix as function of 
m and 4> is now readily available by substituting Eqs. ^ 
and ^ into Eq. ^ . There is an explicit coupling between 
the two order parameters, m and (t> originating from the 
free energy /st — /st(w, 0). Note that the first term in 
Eq. ^ and Eq. ^ lead to a coupHng term which scales 
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as (/)m^. Another possible coupling term is (f)^m that sim- 
ply renormalizes the interaction parameter J of /mix and 
is less interesting. The two-phase coexistence region in the 
(T, (j)) plane is calculated by using the common tangent 
construction to account for the constraint value of the rel- 
ative concentration (/), after minimizing /tot with respect 
to mo and rni. 

By properly choosing the model parameters, we at- 
tempt to reproduce the binary phase diagram of several 
lipid mixtures, the first being that of DPPC/DPPE [18j 
shown in Fig. [21 Since a pure DPPE membrane does not 
show experimentally the rippled phase, its elastic constant 
is chosen as Ca — —2.0 x 10~^^ J. Namely, below the triple 
point Ctr (Fig. [1]) . The DPPE reference temperature 
is calculated from Eq. ([3]) with the same values for the 
Landau coefficients as before, by fitting its main transi- 
tion temperature to the observed one of Tm = 60 °C. On 
the other hand, the elastic constant of DPPC is chosen as 
Cb = —5.2 X 10~^^ J (above Ctr) so that the Ppi phase 
appears between the Lpi and La phases. The reference 
temperature of DPPC is set as = — 15 °C, correspond- 
ing, from Eq. ([3]), to a Lf}' La transition temperature 
T,„ — 36 °C. This temperature is not available from exper- 
iments but Hes in-between the experimentally observed 
Lpi — » Ppi phase transition temperature at 32 °C and the 
Pp' La one at 40 °C. 

All the above chosen parameter values are used to cal- 
culate the binary phase diagram fitting the DPPC/DPPE 
mixture. The relative composition 4> is chosen such that 
^ = corresponds to pure DPPC, while = 1 to pure 
DPPE. In addition, the interaction parameter is chosen 
as J = 1.45 X 10~^° J. When the temperature is relatively 
low, there is a coexistence region between two Lpi phases 

(denoted by L^^} -\-L^p}). This low temperature coexistence 
terminates at a critical temperature, above which there is 
only one Lp' phase for the entire range of composition (f). 
It should be noted, however, that the appearance of the 
critical point between the two "soHd" phases is an arti- 
fact of the model. As the temperature is further increased, 
the Ppi phase appears close to the pure DPPC axis. It is 
bound by a coexistence region with the Lp' phases as well 
as with another (much smaller in its extent) coexistence 
region with the La phase. This latter region is very small 
on the scale of the figure and collapses almost into a Hue. 
These two coexistence regions terminate at a triple point 
temperature, T = Ttr, at which all the three phases co- 
exist. For even higher temperature, T > Ttr, the rippled 
phase disappears and the only coexistence region is be- 
tween the La and Lpi phases. 

The calculated phase diagram is in quantitative agree- 
ment with the experimental one for DPPC/DPPE mix- 
tures (e.g.. Fig. 7 of Ref. [18]). Even the coexistence be- 
tween the two solid phases is alluded in Ref. [18]. A sim- 
ilar type of phase diagram was obtained for mixtures of 
DMPC/DPPS and DPPC/DPPS [13 . In these lipid mix- 
tures, both DMPC and DPPC exhibit the rippled phase, 
while pure DPPS undergoes a direct transition from the 
Lpi to La phases similar to DPPE. 
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Fig. 4. Mean-field phase diagram of a binary mixture of two 
lipids as a function of the relative composition and tem- 
perature T . All parameters are the same as in Fig. [2] except 
Ca = -5.5 X 10"" J and Cb = -6.5 x 10"" J. The lipid inter- 
action parameter is set to be J = 1.45 x 10~^° J. The critical 
point is located at 21.3 °C as in Fig. [2] 
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An even stronger lipid-lipid segregation can be mod- 
eled by increasing the value of the lipid interaction pa- 
rameter to J = 1.85 X 10^^° J, while keeping all other 
parameter values unchanged. In the resulting phase dia- 
gram shown in Fig. [3l the lower coexistence region be- 
tween the two Lfji phases penetrates the other coexis- 
tence regions. Such a phase diagram resembles that of 
DEPC/DPPE lipid mixtures [13j. This combination of 
lipids have a strong segregation tendency because both 
head and tail moieties are different. In the calculated phase 
diagram, there are three triple points indicated by the hor- 
izontal dashed Hues. The rippled phase can now coexist 
with either of the Lpi (regions Pp' -\- L^V or Pp' + L^^,' 

in the figure) as well as with the La phase {Ppi + La^). 
There is also a two-phase coexistence region between the 
two liquid phases (La'' -I- l'^o'). The three-phase coexis- 
tence between the two La phases and the Lpi, and the 
existence of the critical point have been indeed observed 
for DEPC/DPPE lipid mixtures [13]. 

In Figs. m and [U we calculate the phase diagrams in sit- 
uations where each of the two pure Hpids exhibit a rippled 
phase between their corresponding La and Lpi phases. 
Hence, we choose the magnitudes of the elastic constant 
to be larger than |Ctr| of Fig.[lJ Ca = -5.5 x IQ-^^ J and 
Cb = —6.5 X 10~^^ J. The choices of the reference tem- 
peratures and Tg, as well as all other parameters, is 
the same as in Figs. [2] and [H The interaction parameter 
is taken as J = 1.45 x 10"^" J in Fig. IH and have a some- 
what higher value, J = 1.85 x 10~^°J, in Fig. [H In the 
phase diagram of Fig.lH one region of the rippled phase ex- 
tends throughout the entire range of the Hpid composition. 
Above and below this region, the rippled phase coexists 
with the Lpi and La phases. At low temperatures, another 
coexistence region can be seen between the two Lpi phases. 
For a larger interaction parameter J (Fig. [5]), these coex- 
istence regions merge to form a large two-phase region, 

and there are four triple points. Here we see a new coexis- 
ts) 

tence between La and L^/ phases which does not exist in 
Fig. m Although such more complex phase diagrams have 
not been yet reported in experiments, we expect that they 
can be found in the future by properly choosing the lipid 
mixtures. One possible choice of lipids that may give rise 
to such a phase behavior is the DPPC/DMPC mixture. 
Both lipids exhibit the rippled phase between their Lpi 
and La phases. 

4 Summary and discussion 

In this paper we proposed a model describing the liquid- 
solid coexistence region in membranes composed of a bi- 
nary lipid mixture. We addressed in particular the possible 
existence of a rippled phase between the liquid and solid 
ones, and its effect on the global phase diagrams. In order 
to incorporate the possibility of a rippled phase, we use 
a model previously proposed for single-component mem- 
branes by Goldstein and Leibler (GL) [T91I2Q] . Their main 
purpose was to present a model describing the lyotropic 
lamellar phases of lipid bilayers as stacks of interacting 



membranes. Here we have shown explicitly that the GL 
lipid stretching energy of an isolated membrane exhibits 
the sequential phase transitions Lpt Ppr — > La by in- 
creasing the temperature. We further extend the single- 
component model and apply it to binary lipid mixtures 
assuming that the elastic parameter C as well as the refer- 
ence temperature T* depend linearly on the relative com- 
position of the two lipids. The calculated phase diagrams 
are in quantitative agreement with the experimental ones 
for specific lipid mixtures. We have also predicted other 
types of phase diagrams, which have yet to be checked 
experimentally for other lipid mixtures. 

Several points merit further discussion. In Sec. [2l we 
choose the relative membrane thickness as the scalar order 
parameter to describe the membrane phase transitions [191 
[20] . In another model proposed by Chen, Lubensky and 
MacKintosh (CLM), a two-dimensional molecular tilt was 
used as a vector order parameter [34ll35] . In their model, 
a coupling mechanism between the membrane curvature 
and the gradient in molecular tilt was considered. As a re- 
sult, various types of rippled phases with different inplane 
symmetries have been predicted. 

In spite of the success in predicting, e.g., the square 
lattice phase that was found in the experiments [40], the 
CLM model suffers from some deficiencies which makes it 
inappropriate to be used in our study. Within the CLM 
model, the Lpi Ppi phase transition cannot be induced 
only by changing the temperature, while keeping all other 
system parameters fixed. In other words, the phase bound- 
aries in the CLM model is parallel to the temperature axis 
(see, for example. Fig. 9 in Ref. ^5]). This is essentially 
a result of the vector nature of the CLM order parameter 
[41] . Unfortunately, such a phase behavior is not in accord 
with the experimental phase diagrams [27,28j, where the 
Lf3> -Pp' phase boundary depends on temperature. This is 
the main reason why we have chosen to use a scalar order 
parameter giving us the phase diagram as in Fig. [TJ 

The second problem is that the CLM model predicts 
a second-order phase transition between the La and Lpi 
phases. However, the main transition is known to be first- 
order according to various experimental investigations [iHl 
[20] . as was previously mentioned in Ref. [40]. A possible 
way out would be to include the 6th order term of the tilt- 
ing vector in the CLM Landau expansion, this lies beyond 
the scope of the present work. 

The two problems mentioned above do not exist in the 
GL model. However, since the choice of a single scalar as 
the order parameter is a major simplification, many de- 
grees of freedom of the lipid molecules are not properly 
taken into account. For example, one cannot distinguish 
between the similar Lpi and Lp phases because the scalar 
order parameter m cannot model the molecular tilt char- 
acteristic of the solid-like Lp phase. With this choice of or- 
der parameter, one can only distinguish between a dilute 
and a condensed phase. Because the membrane thickness 
does not account for the different symmetries of the La, 
Lpi, and Pp' phases, m cannot represent any crystalline 
order, and the solid phase is not properly described here. 
Hence the appearance of the critical points between the 
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two Lp> phases in Figs. [2] and [4] is an artifact of our model. 
Although this problem is remedied in the CLM model, 
the GL model captures the essential features of the struc- 
tural transitions in Hpid membranes in the presence of the 
rippled phase. Using the latter model is sufficient to re- 
produce even quantitatively some of the global features 
of phase diagrams of binary lipid mixtures, in agreement 
with experiments. 

In Eq. ^ we assumed for simplicity that the rippled 
phase is modulated spatially only in the x-direction. In 
general, one can consider the full (x, y) inplane modula- 
tions as treated in the CLM model [34,35j. This can be 
accomplished by including in Eq. ^ additional gradient 
terms in the y-direction. The resulting phase diagram will 
be more complex because there can be more than one 
rippled phase. Although a two-dimensional square lattice 
phase has been identified for a single-component DTPC 
Hpid [4^, it has not yet been found in binary Hpid mix- 
tures. We also note that the one-dimensional modulation 
of the order parameter expressed by Eq. (|4]) is similar to 

(2) 

the Pg, phase in CLM, with the main difference that we 
include a non-zero average term toq. The physical inter- 
pretation of the actual order parameter in the two models 
differs as is explained above. 

Finally, we discuss why the elastic coefficient C can 
take negative values in Eq. ([2]) as is required to get the 
P/3' phase at equilibrium. In Sec. [21 we have mentioned 
that the coupHng between the conformation of the chains 
and the curvature of the membrane may be a possible 
reason [19]. In general, if there exists a coupling between 
an additional elastic degree of freedom and the gradient 
of the order parameter, the elastic constant C will be re- 
duced. This situation can be expressed by including the 
following two additional terms in the free energy: 

1 „ 9 / dm\ , , 

/e.= 2^P^-7p(^). (9) 

Here p represents some elastic degree of freedom, B is 
another positive elastic constant, and 7 is a coupling con- 
stant. For example, p can be taken as the curvature of the 
lipid/water interface, and B is the bending rigidity of the 
membrane. In the second term, p is linearly coupled to 
the gradient of m which induces a spontaneous value of p. 
Minimizing /d of Eq. ^ with respect to p, we get 

This means that the equilibrium value of p is proportional 
to the gradient of m. By substituting this equation back 
into Eq. the effective elastic constant becomes 

C^C-j^B, (11) 

in Eq. Hence C is reduced by 7^/5 and can reach neg- 
ative values when the coupling constant 7 is large enough 
(or when the elastic constant B is small enough). Although 
it goes beyond the scope of the present paper to suggest 
specific molecular mechanism to account for this added 



elastic degree of freedom, such a mechanism would, in 
principle, explain the negative value of the effective elastic 
constant, C. 
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